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Abstract We explore the predictions of a class of dark 
energy models, quinstant dark energy, concerning the 
structure formation in the Universe, both in the linear 
and non-linear regimes. Quinstant dark energy is con- 
sidered to be formed by gumtessence and a negative 
cosmological constant. We conclude that these models 
give good predictions for structure formation in the lin- 
ear regime, but fail to do so in the non-linear one, for 
redshifts larger than one. 
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1 Introduction 



In a fo rmer publication (jCardone. Cardenas, and Levva 



( 2008^ 1 ) we explored the potential degeneracy existing 



between some dark energy models and f{R) modi- 
fied gravity theories. This means that in principle 
we might reproduce the dynamics of the expansion 
of the Universe with a dark energy model, and then 
obtain an equivalent f{R) model with the same dy- 
namics. In that paper we worked with a class of dark 
energy models which could be dubbed quinstant dark 
energy, as this ingredient is considered to be a com- 
posite of a gwmtessence scalar field and a (negative) 



to ( Cardenas 


et al. 


2002). Cardenas et al. (2002). 


Cardenas et al. 


(2003 


. Cardone. Cardenas, and Levva 


(20081)') for details on quinstant dark energv. 
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In that work we also wrote on the potential of struc- 
ture formation to break away the above mentioned de- 
generacy: the same expansion history might not imply 
the same formation history. Therefore, in this paper 
we give a step forwards concerning the predictions of 
quinstant dark energy for large scale structure forma- 
tion, leaving for the future the exploration of the (po- 
tentially more difficult) predictions of f{R) modified 
gravity theories. 



2 Some characteristics of quinstant dark 

energy 

We investigate spatia lly flat, homogeneous and isotropic 



cosmo logical models ([Cardenas et al\ (|2003[ ) ; ICardone. Cardenas, and I 
( 20081 )) filled with three non-interacting components: 
pressureless matter(dust), a scalar field (/) and a neg- 
ative cosmological con stant A. We first consider th e 
potential introduced bv lRubano and Scudellarol ( 20021 ): 
[Cardenas et al. (2002,): 



V{(j)) oc exp{- 



(1) 



For this potential the following substitutions inspired 
in the Noether symmetric approach allow analytical so- 
lutions: 



uv 



(2) 
(3) 



where a is the scale factor, ui and U2 are functions of 
time. This makes it possible to integrate the Friedmann 
equations exactly. Setting at present (r = 0): a(0) = 1, 
d(0) = 1 and H{0) = 1 we obtained the evolution of the 
scale factor as a function of time: 



_2tjV^ 



i(ai + bitan ujT){a2 + b2tan lut) 



(4) 



2 



with: 



01 = 1 

bi = Ui/U2 

02 = w{2y^-^A - UiU2t) 

62 = 2v/-fiA(3 - U1UJ/U2) + W1U2 + 

w = 2A/-riA/2 



(5) 




(6) 


Wo 


(7) 




(8) 


ax 


(9) 


T 



Fig. [T] shows the parameterisation. The best fit is 
obtained with the foUowing values of the parameters: 



= 0.97944 ±0.00097 
= 0.49939 ±0.00012 
= 0.10254 ± 0.00021 



(11) 
(12) 

(13) 
(14) 



3 Parameterisation of the equation of state 

Because of the rapid evolution of the equation of state 
in our solutions, the conventio n al parameterisations 



This approach allows integration over w to be done 
analytically: 



/(«) 



exp 



(15) 



Huterer and Turneil (ll999l.l200lD:IWeller and Albrecht 



■20011 ): IChevallier and Polarskil (|200ll ): iLinded (|2004 
do not work properly. In order to parameterise 
it we use a 4-parain e ters e quation of state proposed by 
Linder and Huteren (j2005). This is similar to the Kink 
aproa, c h proposed by iBassett. Corasaniti. and Kunz 



and so the Hubble parameter can be written explic- 
itly: 



( 2004|) : ICorasaniti and Copelandl (|2003l) : 
2OO3I) . 



n,noa-^ + (1 - n™o)a-3(i+"'") 

-3tAi 



Corasaniti et aL^ 



1 



-l/r 



3tAw 

1 + (—)!/-) (16) 
ax 



w{a) = wq + 



Wo 



(10) 



Here ut is the value of the scale factor at a transition 
point between w — Wm, the value at the matter domi- 
nated era, and w — wq, the value today. The rapidity 
of the evolution is measured by t 



4 Linear growth of fluctuations 

The dynamics of these models has been successfully 
probed using the standard test of the acoustic peak, 
the dimensionless coordinate distance y(a) and others. 
In this work we must study the role of dark energy on 
the evolution of density perturbations. 
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Fig. 1 Behaviour of the equation of state (here called w) 
versus the scale factor. The red line corresponds to the 
4— parameters approach. The black dotted lines are the 
limit cases, from top to bottom: a) w = 1 stiff matter, 
b) w = dust, c) w = —1 cosmological constant. 



Fig. 2 The growth history(D(a) — S/a) is shown for the 
composite dark energy model(red line). The black line cor- 
responds to the ACDM model 



The perturbation equation reads as: 
S + 2HS ~ AttGpmS = 



(17) 
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It is convenient to study the growth evolution in 
terms of the scale factor a or the characteristic scale 
In a rather than time t. Because in the matter dom- 
ination epoch the solution is (5(a) ~ a, it is useful in 
studying the influence of dark energy to divide out this 
behavior and shift to the growth variable D{a) = |. So 
we can rewrite equation p7p as: 



D" + D' 



5 {In E^y 
a ^ 2 



D 

a 



1 - 



n 



M 



{In E^y 



= 



(18) 



where the prime denotes derivation with respect to 
the scale factor a and E^ = {jj^Y- This equation 
is solved numerically using the boundary conditions 
D{aLSs) = 1 and D'{aLSs) = 0. Here qlss = j+j^ 
denotes the scale factor evaluated at the redshift of the 
last scattering surface. In our calculation we used the 
approximated value of zlss ~ 1089 

The growth history is shown in Fig. [21 more informa- 
tion is given by Fig. [31 which represents the percentage 
deviation {AD {a) = l — D{a)/DA{a)) of the growth fac- 
tor for the composite model of dark energy with respect 
to that for the concordance KCDM one as a function 
of the scale factor. Current deviations are less than 8%. 



Q 




Fig. 3 Percent deviation of the growth variable D from 
the ACDM predictions 

In the linear perturbation theory the growth index 
can be measured by the peculiar velocity field v. This 
quantity is defined as: 



/ = 



din (5 

din a 



(19) 



Nevertheless, it is better to solve a direct equation 
for /, to avoid propagating numerical errors from the 



calculation of 5. For this purpose, one should simply 
use the definition of / across (|20l) and Eq.([T8|) to get 
the evolution equation of the growth index versus scale 
factor: 



/' + 



2 {InE^y 
a ^ 2 



f- 



3n 



M 



2£;2a4 



= 



(20) 
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Fig. 4 The growth index as a function of the scale factor 
for the composite dark energy model(red line). The black 
line corresponds to the ACDM model 




Fig. 5 Percent deviation of the growth index / from the 
ACDM predictions 

This equation is solved numerically using the ini- 
tial condition f{aLSs) = 1 and Eq. p^ . The result is 
shown in Fig. [H To analyse the behavior of / predicted 
by the concordance scenario and the composite dark 
energy model we compute the percentage deviation be- 
tween their predictions as A/(z) = 1— /(z)//a(z). This 
is shown in Fig[5l It is worth noting that Af is always 
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negative, i.e., the growth index of the composite dark 
energy model is larger than the ACDM. Despite this 
result the predicted value of the growth index at the 
survey effective depth z = 0.15 is in very g ood agree 



V 



1/2 r 



X.t: 



(26) 



ment with the estimated value reported in ISen et 
(I2OO6I) 



where we have considered that the spherical overdensi- 
ties evolve in the presence of homogeneous dark energy 
(generally dark energy does not cluster on scales smaller 



than 100 Mpc, see Dave. Caldwell, and Steinhardtl 



5 Non linear evolution of a spherical 
fluctuations 



and also the mass of the forming cluster is con- 
served: 



Pcluster'^ — P cluster ^ta'^ ta 



and 



Another powerful tool to understanding how a small 
spherical patch of homogeneous overdensity forms a 
bound system via gra vitational instab il ity is the spheri- 
cal collapse modelfsee Basilako j ( 20031), Mot a and van de Bm cl^ p, 
j2004VHo rdlou and Bergd |2005l ). lMaor and Lahavl (|2005h .i ~ \pbackg 
IWang (.20061 )). The equation of motion of a spherical 
shell in the presence of dark energy is 



^cluster 



(27) 



(28) 



47rG 



(a 



'cluster 



■ PX.eff) 



(21) 



where pduster is the time varying density inside the 
forming cluster and 



PX.eff = Px + 3px = (3w(r) + l)px 



(22) 



is the effective energy density of the X— component. For 
a cosmological constant case px,eff is constant. Note 
that for w{r) ^ —1 and w{r) ^ —1/3 the density of the 
X— component inside the overdensity patch depends on 
the evolution of the background. If the initial overden- 
sity exceeds a critical value, the overdensity region will 
break away from the general expansion and go through 
an expansion up to a maximum radius (turnaround 
point), then it will begin to collapse until the sphere 
virializes forming a bound system. 

5.1 Rescaled equations and evolution until turnaround 

Let us perform the following transformations: 



a 

ata 



and 



y 



rta 



(23) 



(24) 



where ata is the scale factor of the background and rta 
is the radius of the overdensity when the perturbation 
reaches the turnaround. In a flat cosmology (fc = 0), 
using equations (|16l2ip . the evolution of above quanti- 
ties is given by: 



X ^ Htan\'a^[n{x)x\ 



-1/2 



(25) 



is the overdensity of the forming cluster at turnaround. 
Its value can be calculated by integrating the above 
differential equations (|25l26p using at the same time the 
boundary conditions {dy/dx)x=i — and yx=o — 0. 
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Fig. 6 Evolution of a perturbation collapsing at the 
present time {zcou = 0) in different cosmologies. The values 
Wo — —0.75 and wi = 0.4 were used in the third dark energy 
parametrization. The x— axis shows the time normalised to 
the turnaround time for the ACDM model. 

The variations of a perturbation of a collapsing now 
in different cosmologies are compared in Fig. [51 The 
composite dark energy model reaches turnaround and 
collapse earlier than the other models, this behavior 
is because of the effect of the negative cosmological 
constant which contributes with a positive pressure in 
the whole evolution. The evolution of the overden- 
sity C with collapse redshift is shown in Fig. [71 All 
quintessence models, including the w = wq + wi{l — a) 
parameterisation, show that perturbations collapsing at 
certain redshift are denser at turnaround, relative to the 
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background, than in the ACDM model. These models 
tend at high redshift to the fiducial value of 5.6 for the 
Einstein-de Sitter universe. Nevertheless, the compos- 
ite dark energy model shows an unusual evolution be- 
cause its relative overdensities grow while the redshift 
tends towards a higher values. 
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Fig. 7 Evolution of a perturbation collapsing at the 
present time {zcou ~ 0) in different cosmologies. The com- 
posite dark energy model, upper left side of this figure, has 
an unusual behavior 

This worrying behavior is more evident in Fig. [HI 
This is a first important departure of quinstant dark en- 
ergy from the mainstream interpretation of astrophys- 
ical observations. 
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Fig. 8 Evolution of a perturbation collapsing at the 
present time {zcoii ~ 0) for the quinstant dark energy model. 
The overdensity grows indefinitely a,s Zcoii tends to high red- 
shift. 



5.2 Density contrast at virilization 

In order to obtain more information about the pre- 
dictions of the composite dark energy model from the 
growth history point of view, we need to compute the 
density contrast of the virilised cluster as a function of 
the collapse redshift: 



Pv 



Pbg.vi 



1 + Zta 
1 + ZcoU 



(29) 



where rj = Pvir/ Pta is the ratio of density of the viri- 
alised sphere to the density at turnaround. This pa- 
rameter is very sensitive to the way you consider the 
virialisation process. The usual method to find r/ is 
combining the virial theoren{3 and conservation of en- 
ergy to obtain a relation between the potential energies 
of the collapsing sphere at turnaround and at viriali- 
sation time. In this process the gravitational potential 
energy of the spherical dark matter overdensity will be 
modified by a new term due to the gravitational effects 
of dar k energy on dark matter (see iHorellou and Berge 
(|2005( ) and references therein). 

When dark energy is dynamical, i.e. quintessence, 
the total energy of the bound system is not con- 
served and will contribute as a non-conservative f orce t o 
the dark matter part icles (see Maor and Lahav ( 20051 ): 
Zeng and Gaol (|2005l )). One way to avoid this issue 
in the context of an h omogeneous d ark energy com- 
ponent was proposed in IWanj ( 20061 ). This approxima- 
tion considered that dark matter particles can reach a 
quasi-equilibrium state in which virial theorem holds 
instantaneously because during the matter domination 
era the effect of dark energy is still small. 

Thus assuming dark matter has reached this quasi- 
equilibrium state, its total energy can be computed by 
the virial theorem : 

4g(l + Zw^^^„W - 2(1 + (1 + 3w^,ta)(?)a; + 1 = (30) 

where q = Px,ta/ Pmc,ta characterises the strength of 
dark energy at turnaround. Fig. [S] shows the evolu- 
tion of q on Zcoii for the composite dark energy model. 
The effect of energy nonconservation is not so worrying 
because q will always be quite small 

The variations of the ratio of the virialisation radius 
to the turnaround radius is shown in Fig. 1101 The 
shape of the composite dark energy evolution is over the 
standard Einstein-de Sitter universe(pink dotted line), 
the value of ?7 > 0.5 is due to the fact that the effective 
repulsive force of dark energy on dark mater particles 



^ both components of the quinstant dark energy model have their 
potential energies proportional to B? 
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Fig. 9 Evolution of tlie parameter q respect to tlie virial- 
isation redsliift Zyir- 



allows them to reach the quasi equilibrium with a larger 
radius. This behavior is less marked in the composite 
dark energy model than the rest of the models because 
the negative cosmological constant contributes to the 
evolution with a positive pressure in contrast to the 
negative one of the quintessence. 
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Fig. 10 Ratio of the virialization radius to the turnaround 
radius rj — r^ir/rta as a function of the z^ir- 

Finally, the density contrast at virialisation is shown 
in Fig. [TT] All the models have overdensities that de- 
crease with the value of the Zcoii but the composite dark 
energy has the opposite behavior, so this result is con- 
sistent with the previous one shown in Fig. [7] and Fig. 

m 



Fig. 11 Overdensity A„ir as a function of the collapse 
redshift for some cosmological models. The composite dark 
energy model shows an unusual evolution. 



6 Cluster Number Counts 

Let us now study the cluster abundance in the Universe 
using the Press- Schechter formalisn]!. Following this 
approach the comoving number density of collapsed ha- 
los with masses in the range M and M + dM at a given 
redshift z can be calculated using the standard equa- 
tion: 



dn 
dM 



2 p„ 



IT M a{M,z) 



dlna(M,z) 

■e '■ 



dM 



(31) 



where pmo is the present average matter density of the 
Universe, Sc{z) is the linear overdensity of a perturba- 
tion collapsing at redshift z and (t(M, z) — D{z)aM 
is the rms density fluctuation in spheres of comoving 
radius R containing the mass M. 

We have used the fit provided by Weinberg and Kamionkowski^ 
( 20031 ) to 5c{z) because this parameter shows a weak 
dependance on the cosmological constant. 



5c{z) = 



_ 3(127r) 
20" 



2/3 



- aLogiaVLmiz)] 



(32) 



0.353w2 + 1.044u)3 + lT28w2 + 0.555?i; + 0.131 (33) 



•^it is possible to us e m ore realistic mass function(see 
ISheth and TormenI |l99§) and ljenkins et al\ l|200ll ')') but the aim 
of this section is to draw a comparison between the general predic- 
tions for the cluster number counts of the composite dark energy 
model and ACDM 
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Fig. 12 Evolution of number counts in mass bins witli 
redsliift for objects witli mass witliin tlie range 10" < 
M/{h-^MQ) < 10". 
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Fig. 14 Evolution of immber counts in mass bins with 
redshift for objects with mass within the range 10^^ < 
M/{h-^MQ) < 10^^ Lines are the same as in Fig. 1121 



We use the fit provided in Viana and Liddle ( 19961 ) 
to obtain the quantity 



Cm = CTs 



M 



(34) 



where as is the variance of the over-density field 
smoothed on a scale of size Rs = 8h~^Mpc, its ac- 
tual value is constraine d by the observations of WMAP 
( Komatsu e^"aZI (|2009l) ) andMg = 6 x lO^'^f^^o^'^MQ. 
The index 7 is a function of the mass scale and the 
shape pa rameter (F = 0.149) of th e matter power spec- 
trum (see Viana and Liddld ( 19961 ) for a detailed deriva- 
tion of these quantities): 



7(M) = (0.3r + 0.2) 







2.92 +i log 









(35) 
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Fig. 13 Evolution of number counts in mass bins with 
redshift for objects with mass within the range 10^'* < 
M/ih-'^Mo) < 10'^ Lines are the same as in Fig. 1121 



In order to compare the different models we use the 
normalisation equation (label mod stands for model): 



ScMjz = 0) 

Sc,Lambda{z — 0) 



CT8,A 



(36) 



The value of g8.A = 0.7 96 is also constrained by WMAP 
( Komatsu et al. ( 20091 )) . In this work the effect of dark 
energy on the number of dark matter halos is analyzed 
by computing the number of halos per unit of redshift, 
in a given mass bin 



dN 

Nun = — I dil 

dz 



Ms. 



dn d'^V 



where the comoving volume is 



dM dzdn 



dM 



dzdfl 



dAiz)\l + zf 

' H{z) 



(37) 



(38) 



and dA{z) is the angular diameter distance. 

The number counts in mass bins, Nun = dN/dz, ob- 
tained from ([37]), are shown in Figures. [ni[l31[lll All 
the models exhibit similar behavior: the more massive 
structures are less abundant and form at later times, as 
it should b e in the hierarchical model of structure for- 
mation (see Liberato and Rosenfeldl ( 2006 ) for a further 



discussion in another dark energy parametrization). 

Concerning the quinstant dark energy model, it is 
capable of reproducing the results of the other models 
in a satisfactory way backwards in time up to redshifts 
a bit larger than z = 1 for the three range of mass 
values. Then, it shows abrupt peaks of structure for- 
mation, in a serious departure of the hierarchical model 
for large scale structure. This seems to be caused by 
the unusual equation of state of quinstant dark energy, 
which behaves as stiff matter for redshifts a bit larger 



8 



than one. This would result in enhanced accretion of 
the forming structures, both because of gravitational 
and viscous forces. 



7 Conclusions 

Our opinion is that any dark energy model which 
shows a stiff-like equation of state in the past, dur- 
ing a long period of time, will predict abrupt peaks 
of structure formation, which would be the result of 
enhanced accretion of the forming structures, both be- 
cause of gravitational and viscous forces. This seems 
to be t he case for all studied models of quinstant dar k 
energv (ICardenas et al. I (l2002l).ICardenas et al\ (I2OO2I) 



(|2onah i. 



Cardenas et al\ (|2003f ). ICardone. Cardenas, and Levva 



There are other cosmological models of com- 
posite dark energy having stiff-matter domination as 
an attractor in the past, but usually they would not be 
global attractors, but local. This is the case of several 
models of quintom dark energy. 

Our general conclusions on structure formation pre- 
dictions of quinstant dark energy (QDE) are: 

• QDE makes reasonable predictions for the formation 
of linear large scale structure of the Universe. 

• It reproduces reasonably well the non-linear struc- 
tures from today up to redshifts a bit larger than 
one 

• It fails to reproduce the perturbations in the non- 
linear regime for redshifts a bit larger than one. 

In the future we would study how the dynami- 
cally equivalent models of f{R) modified gravity would 
behave concerning structure formation, and then we 
would have a better understanding on how these obser- 
vations would crack the degeneracy dark energy-/ (i?) 
modified gravity. 
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